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Abstract 

The age of globular clusters inferred from observations depends sensitively on assump- 
tions such as the initial helium abundance and the mass loss rate. A high helium abundance 
(e.g., Y «0.28), as well as an inclusion of helium diffusion and oxygen-enhancement in stellar 
models, can lower the current age estimate for metal-poor globular clusters from Mil. 5 Gyr 
to about 11±1 Gyr, significantly relaxing the constraints on the Hubble constant, allowing 
values as high as 60km/sec/Mpc for a universe with the critical density and 90km/sec/Mpc 
for a baryon-only universe. The uncertainties of a high helium abundance and an instability 
strip induced mass loss near the turn-off in globular clusters are discussed. Ages lower than 
10 Gyr are not possible even with the operation of both of these mechanisms unless the 
initial helium abundance in globular clusters is > 0.28, which can hardly be accomodated 
by indirect inferences of helium abundances in globular clusters. 
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I. INTRODUCTION 



The latest estimate of the age of globular clusters (GCs) is 14.1±1.5 Gyr averaging over 
24 globular clusters JI|, or 15±1.5 Gyr for oldest GCs 0. This sets a lower bound on the 
age of the universe which severely constrains Hq, the Hubble constant, for a A (the vacuum 
energy density) =0 Friedman- Robertson- Walker cosmology For example, if the age of the 
universe is 15 Gyr, for Q (the ratio of the density of the universe to the critical density) = 1, 
the Hubble constant Hq must be less than 45km/sec/Mpc ffl; while for a baryonic universe 
where 0.01//i 2 < Vl < 0.015//i 2 (with h = i^ /100km/sec/Mpc) as constrained by Big Bang 
Nucleosynthesis fH,f§], H < 70km/sec/Mpc. Currently dynamic measurements of the density 
of the universe favor Q > 0.3 f|, which yields Hq < 55 km/sec/Mpc. Direct measurements 
of the Hubble constant, however, yield values that range widely from 40km/sec/Mpc to 100 
km/sec/Mpc [|]. 

A reconcilation of the direct measurements of the Hubble constant and that inferred from 
the age of the globular clusters may well require not only improvements of the distance scale 
systematics but also continued investigation of the implicit assumptions in the inferred age 
of globular clusters. Otherwise, the conflict may prove to be an insurmountable problem for 
Q — 1, A = cosmology models. 

It is noted that the estimate of the age depends sensitively on many uncertain aspects 
of globular clusters, such as the actual turn-off luminosity (distance) determination, the 
metallicity, the oxygen to iron ratio, helium diffusion in stars, the initial helium abundance, 
and an instability strip induced mass loss inferred from lithium observations of stars [0. 
The uncertainties of the turn-off luminosity, the metallicity and the oxygen to iron ratio 
have been widely discussed in the literature PHIT|, and has been included in the current age 
estimate of 14.1 ± 1.5 Gyr JT|. The helium diffusion in stars is also estimated to lower the 
age of GCs by 5%-15% [|T^-|T4||. Uncertainties in the initial helium abundance and the mass 
loss rate of stars near the turn-off point have been explored previously [|7|,§ . But since these 
systematic uncertainties should not be added in quadrature like statistical uncertainties, a 
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full exploration including uncertainties in the helium abundance, mass loss, helium diffusion 
and oxygen-enhancement is essential to assess their influence on the age estimate of globular 
clusters. 

The aim of this paper is to explore the uncertainties in the helium abundance and 
mass loss in globular clusters both theoretically and observationally, and calculate the age 
estimates for globular clusters under these uncertainties as well as oxygen-enhancement, 
using stellar models with helium diffusion. Section 2 discusses the methods used in dating 
the age of globular clusters and their dependence on various parameters. In this section, we 
also present our formula that relates the age of globular clusters to their parameters, based 
on our stellar models with helium diffusion. Section 3 discusses uncertainties in the current 
inferences of the helium abundance in globular clusters due to uncertainties in both stellar 
evolution theory and stellar observations. In particular, it argues that the current inferences 
may allow a helium abundance as high as 0.28 in globular clusters, which is significantly 
higher than the primordial value of 0.24. Section 4 provides possible scenarios that can 
enhance the helium abundance in globular clusters from the primordial value to a higher 
one. Section 5 then discuss a possible instability strip induced mass loss near the turn- 
off region of globular clusters, and its effect on the luminosity functions (LFs) of globular 
clusters and their age. The observed luminosity function of globular clusters thus constrains 
the size of the mass loss effect. The effect of a mass loss combined with a high helium 
abundance is also investigated. Section 6 summarizes the previous sections and concludes 
that by considering the uncertainties in the helium abundance and mass loss in globular 
clusters, as well as helium diffusion and oxygen enhancement, the age of metal-poor clusters 
can be as low as 11 ± 1 Gyr. But ages lower than 10 Gyr for metal-poor clusters cannot be 
consistent with observed properties of globular clusters. 
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II. AGE DATING OF GLOBULAR CLUSTERS 



Two methods are commonly used to estimate the age of a globular cluster (see review 



of ref. [15]). One is to fit theoretical isochrones (temperature- luminosity curves for stars 
of a given age) to the shape of a globular cluster on the observed color-magnitude dia- 
gram, using a distance modulus determined from fitting either a zero age main sequence 
(ZAMS) or a zero age horizontal branch (ZAHB). In the fitting, it is necessary to transform 
a temperature-luminosity diagram of theoretical isochrones into a color-magnitude diagram 
with an interstellar reddening correction [IT|. Besides the uncertainties mentioned in the 
introduction, this method is very sensitive to the transformation between the luminosity- 
temperature diagram and the color-magnitude diagram, the fiducial sequence of globular 
clusters extracted from their scattered color-magnitude diagram, the subjectivity in fitting 
isochrones to the fiducial sequence, and the effectiveness of a constant mixing length in de- 
scribing the convection. It also depends sensitively on the adopted interstellar reddening if 
using ZAMS fitting to obtain distances. 

The difficulty in improving the fitting method lies mainly in the fact that the color of 
both observed GCs and theoretical isochrones on the color-magnitude diagram is hard to 
determine precisely. The color of observed GCs depends on the interstellar reddening which 
can be very uncertain; while the color of isochrones is sensitive to the treatment of convection 
in stellar models. Figure 1 shows isochrones with similar parameters but different mixing 
lengths (parametrized by a, the ratio of the mixing length to the pressure scale height) in 
the convection zone. Both isochrones have a similar turn-off luminosity, but the color shift 
due to a different choice of a is manifest, especially on the upper subgiant and the red giant 
branch of isochrones. Although the uncertainty in choosing different mixing lengths may be 



partially resolved by calibrations of ZAMS stars [|lq| , it still persists since mixing lengths 
may change in different evolution phases, such as the subgiant branch and the red giant 
branch. 

A method that is not sensitive to the color determination is to relate the age of a GC 
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with the difference in luminosity between the GC's turn-off point and horizontal branch, 
regardless of its detailed morphology. Such a relation has been shown by Iben and Renzini 
HI to be 



log(^) = 8.497 - 1.88(Y - 0.24) - 1.44(Y H b -Y)- 0.088 log Z + 0.41 (M^ - M ] 

(1) 



RR\ 
bol )i 



where Y is the initial helium abundance of GCs, Z is their metallicity, Yhb is the helium 
abundance at the envelope of the HB stars, and M^J are the bolometric magnitudes 

of the turn-off point and the RR Lyrae stars respectively. The coefficients in the equation 
have been refined over the years, but the resultant age hasn't changed by more than 25%. 



For example, a similar relation based on Sandage's calibration of RR Lyrae stars [18] and 



the oxygen-enhanced isochrones of Bergbush and VandenBerg |TlJ, assuming Y « 0.24, is 

m 



log(y^) = 8.773 + 0.017[Fe/H] + 0.39(M™ - M* R ), (2) 



where M v denotes the visual magnitudes. Equation (Q) tends to yield an age that is ~ 20% 
lower than eq. (|I|) for metal-poor clusters. The method of relating the age of GCs and their 
turn-off luminosities is independent of the interstellar reddening, the subjectivity of fitting 
isochrones to the observed GC sequence and the mixing length adopted in stellar modeling. 
It still suffers from the uncertainties mentioned in the introduction. 

In order to test the sensitivity of GC ages to the helium abundance for stellar models 
with helium diffusion, we constructed a series of models between 0.60M Q and 0.90M Q , with 
an interval of O.O5M . Their metallicities are chosen to be Z=0.0002 (10~ 2 Z & ), 0.0006 and 
0.002. The helium abundance of the models is set to be 0.24 and 0.28. We adopted a mixing 
length of a = 1.69. All models were started as pre-main sequence objects on the Hayashi 
track, and evolved to helium flush, or 25 Gyr old (which is enough for GCs around 15 Gyr 
old), whichever comes earlier. We used Livermore opacity tables for temperatures above 
6000K |L9]], and Los Alamos opacities below 6000K [pOR . The code we used is based on 



Dearborn's stellar code |2~Ifl, with an inclusion of helium diffusion using the prescription of 



Bahcall and Loeb 1221. From these models, we found that 



log(-j^) = 8.11 -0.12 log Z + 0.40JM££ -0.80(y- 0.24). (3) 



Our models agree fairly well with previous helium diffusion calculations [I5-H|. 

Since the helium diffusion affects little the evolution of HB stars due to their short life 
span, we can just take the result of HB calculations without including helium diffusion 

HIH: 



M RR (bol) = 0.943 - 3.5(F H b - 0.30) + 0.183 log Z. (4) 
Combining eq. (§) and (f|), we get 
log(^) = 8.57 - 0.047 log Z - 2.20(F - 0.24) - 1.40(Yhb -Y) + 0.40(M^ ) - M ] 



RR\ 
bol I 



(5) 

This result compared with non-diffusion calculations based on the same code [|7| lowers the 
age of metal-poor clusters by about 8%. It yields ages that are generally 15% to 25% lower 
than eq. ([]]), depending on Z. 

We couldn't obtain oxygen-enhanced models with Livermore opacity tables. But it has 
been shown that for Los Alamos opacity tables, oxygen-enhanced models with a metallicity 
Z (deduced from [Fe/H] only) can be approximated by standard solar-scaled abundance 
models with a metallicity Z ef[ , provided that |25|] 



Z cS = (0.64/ + 0.36)Z, (6) 

where / is the enhancement factor of oxygen. If this also holds true for Livermore Opacity 
Tables, oxygen-enhanced models will yield the same equation as eq. (H), as long as Z e s is 
substituted for Z. Commonly / is taken to be 3 to 5 in metal poor clusters [|1(]], therefore 



Z c s 3Z, which will further lower the age estimate from eq. (BJ) by about 5%. 
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Yhb—Y is estimated to be 0.01-0.02 in models without helium diffusion ||17|, and —0.005- 
0.005 in helium diffusion models Therefore a strong dependence of the age on 

the helium abundance Y is manifest in both eq. ([!]) and eq. (|5|). 

It is conventionally assumed that Y ~ 0.24 in GCs [fl],[n]|, in agreement with the ob- 
served primordial helium abundance Y = 0.235 ± 0.01 [27] and the prediction of Big Bang 
Nucleosynthesis [HQ. But both eq. fll]) and eq. (|J) show that a higher Y (but < 0.28) and 
can significantly reduce the age estimate of GCs. For example, a Y of 0.28 will reduce the 
age by as much as 20% with respect to a Y of 0.24, according to eq. (|l|) and (||). Even a Y 
of 0.26 will lower the current age estimate by 10%. 

Determining the age of GCs with the isochrone fitting method has a similar Y sensitivity 
as in eq. ([!]) or fl5|) if a distance modulus from ZAHB fitting is used. If a distance from ZAMS 
fitting is used in the isochrone fitting, the effect of a higher Y in GCs on their age is not easy 
to quantify. Figure 2 shows the comparison of a Y = 0.28 isochrone and Y = 0.24 isochrones 
in fitting the fiducial sequence of NGC6397 (adopted from PB"|). Both isochrones include 
helium diffusion but don't include oxygen-enhancement due to the adoption of Livermore 
opacities in most of the temperature range, and have a. Z = 0.0002. NGC6397 is believed 
to have [Fe/H]=-1.91 ||T|. The reddening and distance modulus adopted for the cluster are 
0.11 and 12.35 for the Y = 0.24 isochrones, 0.14 and 12.45 for the Y = 0.28 isochrone. The 



distance modulus estimated from calibrations of RR Lyrae stars is about 12.5 p8| . Fig. 2 
shows that while NGC6397 has an age of ~ 15 Gyr when fit with the Y = 0.24 isochrones, 
it is only ~ 12 Gyr old when fit with the Y = 0.28 isochrones. This is in general agreement 
with eq. (Bh from which a Y of 0.28 yields an age 20% younger than a Y = 0.24. 



III. THE MEASUREMENT OF HELIUM ABUNDANCES IN GCS 

The surface temperatures of stars in globular clusters are too cool to have detectable 
helium lines. Therefore, the Y in GCs cannot be measured directly from spectra of stars. 
Instead, Y is indirectly determined from comparisons between certain theoretically Y de- 
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pendent parameters of GCs deduced from stellar models and the measurement of such pa- 
rameters from GC observations. 

One such parameter is the relation between the effective temperature and period of RR 
Lyrae stars at the blue edge of the instability strip |29j . Comparison of this relation between 
canonical stellar models (i.e., no rotation, no diffusion, no mass loss, etc) and observations 
yields Y HB = 0.27±0.04 @. A similar comparison is done to the A parameter, which relates 



the effective temperatures and periods of general RR Lyrae stars [j3l|. This comparison 



usually yields a high Yhb of 0.26-0.33 [[3TJ, and leads to an anti-correlation between the 



metallicity of a cluster and its helium abundance, which can hardly be explained with any 



chemical evolution theories |31|j32| . This anti-correlation problem may vanish if one assumes 



a varying core mass in HB stars due to non-canonical effects (for example, rotation 



and/or a varying CNO abundance | 3~1|| . Both of these comparisons, however, involve the 
determination of the effective temperatures of RR Lyrae stars, which, as mentioned above, 
cannot be confidently inferred from their observed colors due to the uncertainty in the 
interstellar reddening and stellar modelling. 

Another parameter is R, which is the ratio of the number of stars on the horizontal 
branch to the number of stars on the red giant branch with luminosities higher than the 
mean luminosity of RR Lyrae stars. Stellar models show that R is very sensitive to Yhb and 



thus Y. Canonical stellar models yield a R-Y relation |24|,|34] 



Y = 0.186 + 0.37 log R, (7) 

assuming Yhb — Y ^ 0.02, a result from canonical models that exclude helium diffusion. 
Eq. (0) yields an average Y of 0.23 ± 0.02 for 15 GCs [24|. However, if other non-canonical 



effects such as rotation operate in GC stars to increase M c , the helium core mass of HB 
stars, estimate of Y can be higher than canonical values. For example, a core rotation 
rate of 2xl0~ 4 rad/s, a typical main sequence value, will increase M c by O.O3M PD| , and 
increase the estimate of Y from the R comparison by 0.01, since [ |3~3[ | 

dY/dM c = 0.4. (8) 
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Assuming a larger and varying M c , Caputo et al. obtained an average Y of 0.24 ±0.01 when 
combining the i?-method and the A-method ||34|| . Furthermore, as we have shown, if we 



include helium diffusion, Yhb — Y will be only about 0.005. In addition, M c will be slightly 
increased by O.OO5M with respect to that of canonical models due to the helium settlement 
into the core, which results in a 0.002 increase in Ihb- The resultant Y will therefore be 
0.02 higher than the value inferred from eq. (0), just from the inclusion of helium diffusion 
in stellar models alone. 

The luminosity spread of the horizontal branch serves as another parameter that can 
be used to determine the helium abundance. A larger Ihb yields a broader luminosity 
spread among the horizontal branch stars, according to HB models |p3[| . Such comparisons 



for 47 Tuc and M15 suggested a Ihb close to 0.24 [g5|,[3(J. But again, possible systematic 
variations in M c unaccounted for in canonical models may affect the result. For example, 
to bring the Y values from the A-method and the i?-method into agreement, it is shown by 
Caputo et al. that M c in M15 has to be increased by 0.04M Q relative to its canonical value 
P4j |. According to the calculations of Sweigart and Gross a 0.04M Q increase in M c will 
narrow the predicted luminosity spread of a horizontal branch appreciably and thus allow a 
higher >hb given an observed luminosity width of HB. 

A A paramter is also proposed to determine Y, where A is defined as the magnitude 
difference between RR Lyrae stars and the main sequence at (B — V)q (the corrected B — V 



color)=0.7 p7 ]. Comparison with the measured A yields a Y of 0.2 to 0.3. This method 



depends on the determination of the metallicity of GCs, the color determination and 
hasn't accounted for effects such as helium diffusion. 

Despite these many methods of determining Y for GCs, a statistically sound upper limit 
on Y is lacking. However, a Y > 0.28 is very hard to reconcile with most of the above 
determinations, especially the comparisons of the R parameter and the luminosity spread 
of HB, even when aforementioned non-canonical effects are taken into account. 

A maximal helium abundance of 0.28 in GCs seems also to be consistent with that 
inferred from luminosities of the HB stars. Ihb is related with the luminosity of RR Lyrae 
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stars by eq. (£|) in canonical models. Current calibrations of the visual magnitude of RR 
Lyrae stars yield a 0.2-0.3 magnitude spread for different methods, and an uncertainty of 
0.2 magnitude within each method P, p^8| , |38| , p9| . which translates into a 0.06 uncertainty in 
Vhb- In other words, a Y of 0.28 can be easily accomodated for individual clusters, given a 
primordial value of 0.24. It should be noted that the observation of field RR Lyrae stars in 
the disk of the galaxy may not serve a valid constraint on Yhb in globular clusters if only 
stars in GCs are enriched in helium. For example, RR Lyrae stars in seven Large Magellanic 
Cloud clusters are on average 0.28 magnitude brighter than field RR Lyrae stars f2|, which 
can be consistent with an enhancement of helium abundances only in globular clusters. 



IV. SCENARIOS TO ENHANCE THE HELIUM ABUNDANCE IN GCS 

A helium abundance higher than the primordial value Y p in GCs is possible to achieve 
under certain scenarios. The simplest scenario is the enrichment by the first generation 
supernovae (SNe) [fF|,|4T)f Measurements of low metallicity HII regions in irregular galaxies 
and extragalactic HII regions showed that their helium abundance Y correlates with their 
oxygen abundance [O/H] by [p7) , [4l[ 



Y ^Y p + 130[O/H] for [O/H] < 2.5 x 10" 4 , (9) 

i.e., dY/dZ w 6. For GCs with [Fe/H] = — 1, their oxygen abundances are about 1/3 of the 
solar value ||, or [O/H] « 3 X 10~ 4 . If GCs follow the same helium-oxygen correlation 
as those low metallicity HII regions, they have Y ~ 0.27, which is significantly above the 
primordial value. Of course, it is debatable whether GCs follow the same K-[0/H] relation 
as low metallicity HII regions, although both of them have similarly low metallicities. The 
sun, on the other hand, has Y « 0.28 and [0/H]= 8 x 10~ 4 , which would require a different 
relation Y m Y p + 60 [O/H] (i.e., dY/dZ = 3). Therefore, solar- type disk stars must have 
had quite a different chemical history from low metallicity gases, like those of GCs and 
HII regions. It is worth noting that VandenBerg, Bolte and Stetson El found a 2-3 Gyr 
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age spread among GCs with [Fe/H]> —1.6 and no significant age spread among GCs with 
[Fe/H]< —1.6. If an enhancement in the helium abundance correlates with an increase 
of metallicity in GCs, as suggested above, the age spread among metal-rich GCs may be 
explained at least in part by the different enrichment of helium in different clusters. 

Exotic models, such as Pop III stars may enrich helium but not heavier elements. 



It is calculated that if the oxygen cores of these Pop III stars are larger than ~ 1OOM , 
the cores will collapse completely into blackholes without ejecting heavy elements. However, 
processed helium in Pop III stars can be deposited into the interstellar medium through mass 
loss and mixing of Pop III stars. The maximal amount of helium enrichment is AY = 0.17 



for an initial helium abundance of 0.24, if all gas is processed once in Pop III stars |f43| , [44 



Therefore, in order to enhance the helium abundance in a 10 6 M G) cloud from 0.24 to about 



'e 
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0.28, at least 2.5 x 10 M gas has to be processed through Pop III stars. This requires quite 
a different mass function from the M a power law with a = —2.35 observed today ||45|| . Since 



the metallicity of the cloud is low, the required mass function could have a > —0.5, so that 
there were no more than ~ 100 SNe in a 10 6 M Q cloud that would deposit metals into the 
cloud. 

An enhancement of the helium abundance may also arise during the formation of globular 
clusters due to the difference in the first ionization energy between the hydrogen and the 
helium (13.6eV vs. 24.6eV). In certain GC formation scenarios, this difference may enable 
helium to recombine while hydrogen are still heavily ionized. As a result, helium may have 
a slower sound speed than hydrogen and is therefore more prone to gravitational collapse. 

One class of models of forming GCs, first proposed by Peebles and Dicke [|!J, argues that 



GCs form from primordial baryon density fluctuations after recombination, which have a 
Jeans mass of a typical GC mass, 1O 6 M . According to the standard recombination picture 
of the universe, hydrogen recombined at T ~ 3000K (redshift z ~ 1000) while helium 
recombined at T ~ 10 4 K or z ~ 3000 (inferred from Fig. 5 of ref. In between z ~ 1000 

and 3000 there existed a period when helium recombined and primarily interacted with 
atoms, ions and electrons, and hydrogen were still mostly ionized and primarily interacted 
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with background photons. Therefore helium can collapse earlier than hydrogen due to its 
lower sound speed. 

In the period after helium recombination and before hydrogen recombination, the helium 
falls quickly into the potential well of dark matter, following the density fluctuation of dark 
matter ||49|| . As hydrogen recombines, the density fluctuations in the hydrogen and the 
helium begin to grow jointly. The overdense region, however, already has an enhanced 
helium abundance from the in-fall of helium early on. Assuming the density fluctuation in 
dark matter at z ~ 1000 to be Sp/p and the primordial helium abundance Y p = 0.24, the 
enhanced helium abundance in the overdense region at z < 1000 is 

y p (l + Sp/p) + 1 - Y p 

To enhance the helium abundance to Y — 0.26 on the 10 6 M Q scale requires Sp/p ~ 0.1 on 
this scale, which means Sp/p on the 1O 6 M scale will become nonlinear at z ~ 100. This 
cannot be realized in the standard cold dark matter (CDM) model, in which structures 
on this scale become non-linear at z < O(10) but may be possible in models where 
primordial seeds such as topological defects provide nonlinear fluctuations on small scales 
PH^2], A Sp/p of 0.1 at z ~ 1000 with ~10 6 M Q scales would correspond to fluctuations 
in the microwave background radiation on scales that are less than 8 arc minutes and thus 
would be smoothed by the finite size of the surface of last scattering. 

In models that generically form ~ 1O 6 M cold clouds in proto-galaxies |56]j57[] , a higher 



helium content may arise if Pop II stars form in the shock waves from the first generation 
supernovae in the cloud (simplified in Figure 3) |58[. The temperature of these clouds can 
stay around 10 4 K for a long time due to a sharp drop in the cooling function of the gas 
cloud at 10 4 K [|B|. If a strong UV photon background exists (which is likely due to the 
presence of supernovae shocks), hydrogen can be partially ionized at this temperature and 
interact primarily with the UV photons in the cloud, but helium remains mostly neutral and 
interacts primarily with atoms, ions and electrons. The sound speed in helium, cn e , is then 
~ 10 km/sec (for a temperature of ~ 10 4 K), while the sound speed in hydrogen, ch, can be 
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much larger. When the velocity of a supernova shock wave drops below ch, the shock no 
longer compresses hydrogen in the cloud into the shock front but continues to do so with 
helium. The enhanced helium abundance is then J7| 

y« (11) 

where M s ^ is the mass swept by the shock, M tot is the total mass of the cloud, Y p is the 
helium abundance before the shock, 5 is the thickness of the shock front and I is the thickness 
of the cloud after being compressed by the shock wave. If the SNe rate is > 10~ 6 yr" 1 per 
cloud, the shock wave can propagate through the entire cloud at a speed larger than Ch c 
|7|j40[|. Then M s h ~ M tot . If I ^> S, i.e., the density inside the shell is much higher than 
the density outside, Y can be significantly larger than Y p . Therefore stars formed from the 
compressed gas at the shock front may have a higher helium abundance than the primordial 
value. 

We have briefly discussed several scenarios to enhance the helium abundance in GCs. 
There may be other possible scenarios, such as a magnetic field, that may also segregate 
ionized gas from neutral gas. If these scenarios operate primarily in GCs, they will not 
have a direct effect on HII regions where the primordial helium abundance is measured, 
nor will they affect the helium abundances in galactic disks. Helium-enriched stars from 
disrupted GCs may only constitute a small population of halo stars. Therefore, there may 
not be evidence in disk stars that support helium enhancement processes. Only more refined 
stellar evolution models and GC observations can better determine the helium abundance 
in GCs. Once again, a spread in Y after these helium enhancement processes may be a 
candidate to account at least in part for the spread in age estimates for different GCs which 
assume a universal Y. 



V. THE UNCERTAINTY IN MASS LOSS 



Possible mass loss near the turn-off region of GCs is another uncertainty that affects the 
age estimate of GCs [|7|||[59l . It is motivated by the closeness of GC turn-off temperatures 
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and the small temperature range 6600±200K, where Pop I F-type stars are observed to have 



severely depleted lithium abundances (the so-called lithium dip) pOySIfl . One explanation 
of the lithium dip is an instability strip induced mass loss ||62|| . 

Lithium burns rapidly at temperatures above 2.5xl0 6 K. It is therefore completely de- 
stroyed at layers with temperatures above 2.5xl0 6 K [p5 |. At the temperature range of F 



stars in which the lithium dip occurs, the base of a star's outer convective zone cannot reach 
this lithium-burning temperature. Therefore, its surface lithium abundance should be its 
primordial value (although some negligible depletion occurs in the pre-main sequence [j34H). 
The mass loss at the surface, however, leads to a deeper penetration of the convective layer 
into the star's radiative interior. When the convective zone reaches down to layers with no 
lithium, the surface lithium abundance can be depleted due to mixing |62] , |65|] . The shape of 
the observed lithium dip on the [Li/H]-Temperature plane can be reproduced quite well in 
the mass loss mechanism with a mass loss rate |62| 



M > 7 x Kr u M yr _1 . (12) 

But the narrowness of the observed Pop I lithium dip constrains any significant mass loss to 
a restricted temperature range between 6500K and 6700K. On the other hand, the absence 
of a severe beryllium depletion in this temperature range in Pop I stars limits any mass loss 
rate to 



M < 1 x KT^Meyr" 1 , (13) 

which coincides with the limit from observations on ionized stellar winds of F type main 
sequence stars p5 |. Alternative explanations of the lithium dip include radiative diffusion 



pTfl , and rotation- induced mixings that leads to a more extensive of mixing between the 
surface lithium and the depleted lithium in the radiative interior |68| . 

The same instability strip which may induce mass loss in Pop I F-type stars may also 
operate in Pop II stars in a similar temperature range, 6500-6700K. These stars generally 
have an undepleted lithium layer of about O.O4M at ZAMS. Given a maximal lifetime of 
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several Gyr they spend in the temperature range, Dearborn et al. have shown that if the 
mass loss rates are > lO _11 M /yr, it should produce a similar lithium dip in Pop II stars 
||69|| . Several halo stars in the vicinity of this temperature range have been observed to 
have depleted lithium, but the existence of a lithium dip in Pop II stars is still far from 
established ||7U|| . It is to be noted, however, that even an absence of the lithium dip in Pop 
II stars doesn't automatically exclude the instability strip induced mass loss, as long as the 
mass loss rates are < 10 _11 M Q /yr. This instability strip induced mass loss does not affect 
the lower main sequence and so avoids the arguments against large-scale mass loss through 
out the entire lower main sequence f6~2fl . 

The position of the instability strip on the Pop II main sequence is quite uncertain. If 
mass loss operates in stars in the Pop I lithium dip, by a simple extrapolation, mass loss may 
also operate in a similar temperature range around 6500-6700 K, which is in the vicinity 
of the turn-off temperature of GCs ||. Observations of low metallicity blue stragglers, on 
the other hand, seems to indicate that the instability strip is bluer than the turn-off point 
of GCs with similar metallicities [7T|. However, blue stragglers may be too pathological to 
draw conclusions on the main sequence of GCs. 

Another justification for instability strip induced mass losses in the turn-off region of 
GCs may come from the red edge of RR Lyrae stars and the calculated dependence of this 
red edge on Y (the envelope helium abundance), Z (the metallicity), M (the mass) and L 
(the luminosity). Observations of RR Lyrae stars found the red edge of the instability strip 
in HB stars lies at logT e ~3.81 for M15 and logT e ~3.79 for M3. Although there hasn't 
been a theoretical calculation of the position of the instability strip for Pop II stars, it has 
been shown for Pop I cepheids |72 
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log Tfre = (const. + 0.106Y) + (0.163 - 0.392F - 0.592Z) log(M/M ) 

+ (-0.074 + 0.074F - 0.437Z) log(L/L ), (14) 

where Tfre is the effective temperature of the foundamental red edge. If we take 
eq. (0), and take Y RB « 0.25, Z < lO" 3 , M/M » 0.65 for RR Lyrae stars, Y < 0.1 
TB| , |26| , M/Mq w 0.75 for main sequence turn-off stars, and L RR /L TO ~ 1.4 we get 



logT FRE (RR) — logT FRE (TO) ~ —0.05. Therefore, the instability strip at the main se- 
quence turn-off luminosity has a red edge of logTp R E ~ 3.85, which is very close to the 
turn-off region of metal-poor globular clusters (T e ~ 3.80 to 3.83 |l3,26|). It is therefore 



conceivable that a consistent instability strip calculation for Pop II stars can extend the 
instability strip to the turn-off region. In this regards, theoretical calculations of the insta- 
bility strip for Pop II stars are necessary to extrapolate the red edge accurately from RR 
Lyrae stars down to the turn-off luminosity. 

When mass loss occurs at the temperature of the turn-off, the small reduction in mass 
causes the model to turn off at a lower luminosity with respect to models without mass 
losses, since the turn-off luminosity of a main sequence Pop II star decreases as the mass 
of the star decreases Stars also appear to spend an increased fraction of their life 

beyond the turn-off (since less massive stars always evolve slower than more massive ones) 
||. The GCs will then (incorrectly) appear older due to their lower turn-off luminosities, 
according to eqs. (§, (|) and (g) ||,||. 

Several authors J?]|| have calculated the effect of a mass loss rate ~ 10~ n M o yr 1 in 
the temperature range of 6500 ± 200K and 6600 ± 200 K on the evolution of GCs. Such 
assumptions result in GCs of 11 Gyr to 13 Gyr old looking 2~ 3 Gyr older. To see the 
effect of the mass loss on our helium diffusion models, we repeated our model calculations 
in section 2 with a mass loss rate 

M= lO- 11 M yr- 1 exp[-(^^) 2 ], (15) 

for both Y = 0.24 and Y = 0.28. Figure 4 shows the resultant turn-off visual magnitude 
vs. age relation for different cases. We find that a mass loss rate of eq. ([15]) will only affect 
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isochrones younger than ~ 12 Gyr for Y = 0.24 and ~ 15 Gyr for Y = 0.28 significantly, 
decreasing their turn-off visual magnitude to the level of about 1 Gyr older isochrones 
without mass loss. Therefore, for a GC age of 14-15 Gyr by conventional estimates, or 
about 11-12 Gyr if assuming Y = 0.28, the decrease in age due to a mass loss rate of 
eq. ( ]T5| ) alone is at most ~ 1 Gyr. This age reduction due to mass loss is smaller than 
the 2-3 Gyr reduction calculated previously from stellar models without helium diffusion 
The reason may lie in that in helium diffusion models, mass loss can strip the outer 
helium-poor layer of stars and thus reduce the effect of helium diffusion on stellar evolution. 
For example, a O.8M \ow-Z model with no mass loss has a surface helium abundance of 
less than 3% at turn-off due to helium diffusion, but a similar model with a mass loss of 
eq. (|15"D has a surface helium abundance of more than 20%. 

Since it has been shown that a high helium abundance of 0.28 can decrease the age 
estimate for GCs from ~ 14 Gyr to about 11 Gyr, from figure 4 it can be seen that a 
combination of Y = 0.28 and a mass loss rate like that in eq. ([T5| ) can further decrease the 
GC age to as low as 10 Gyr. However, to what extent the mass loss mechanism can lower 
the age estimate of GCs is not only determined from its turn-off magnitude-age relation but 
also constrained by the calculated GC luminosity function under mass loss compared with 
observations ||8|,|65||. 

It has been pointed out that mass loss near the GC turn-off point can significantly affect 
the luminosity function of GCs, since the longer life time that stars spend after their turn-off 
due to mass loss results in a bump on the LF in the subgiant region |8|j65[|. In figure 5(a) 
and 5(b), we plot the predicted luminosity functions of GCs with a mass loss rate of eq. ([15]) 
(assuming a Salpeter initial mass function |45[]) for both Y = 0.24 and Y = 0.28, compared 



with observations compiled by Stetson |73" |. For Y = 0.24, figure 5(a) shows that a 12 Gyr 
mass loss LF can be ruled out by the absence of any large size bump at the subgiant region 
in the observation, but a 13 Gyr mass loss LF may still be allowed. For the Y = 0.28 case, 
figure 5(b) shows that the 11 Gyr mass loss LF can be confidently ruled out, but the 12 
Gyr mass loss LF may not be ruled out easily. Therefore, due to the constraint from the 
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luminosity function of GCs, mass loss essentially cannot lower the age of GCs to below 12 
Gyr even when combined with a high helium abundance. This result should also hold if we 
consider oxygen-enhancement in GCs, or if we use another power-law initial mass function, 
since it has been shown that the LFs of GCs in the turn-off and subgiant region are not 
sensitive to both of the two factors J75[ . 

It is interesting to note from figure 5(a) and (b) that mass loss can improve the fit at 
the red giant region of LF between models and observations. But it may be premature to 
draw any conclusion from comparison to this uncertain portion of the data. 

Isochrone fittings also seem to indicate that even with a combination of high Y and mass 
loss, an age lower than 12 Gyr is hard to reach. Figure 6 shows Z = 0.0002 isochrones with 
eq. fljjD and Y = 0.28 at 12 and 13 Gyr old, overlaid on the fiducial sequence of NGC6397, 
as well as a standard isochrone without mass loss and with Y = 0.24 at 15 Gyr (the same 
one as in Fig. 2). Clearly, an isochrone in between 12 Gyr old and 13 Gyr old with Y = 0.28 
and mass loss may fit the observation acceptably (while the 12 Gyr isochrone with Y = 0.28 
and no mass loss fits fairly well in figure 2). Therefore, effects of a higher Y and a mass loss 
on GC ages do not add up linearly in the isochrone fitting method, and an GC age lower 
than 12 Gyr seems still unachievable. 

Applying our mass loss calculations to Pop II stars, we find that with a peak mass 
loss rate of 1O -11 M yr l , severe lithium depletions occur at stars between ~ O.8M and 
~ I.OMq (similar to the conclusion of Dearborn et al. ||69|| ). Therefore, in isochrones older 
than ~ 14 Gyr for Y = 0.24 or older than ~ 12 Gyr for Y = 0.28, no stars with lithium 
severely depleted from mass loss will be observed since they already evolved away from the 
red giant branch. That is, if further observations find no evidence for a mass loss induced 
lithium depleted in Pop II stars, it may simply indicate that our universe is too old for 
those lithium depleted Pop II stars (if there were any) to survive. In isochrones younger 
than these ages, lithium-depleted stars should be observed in either their main sequence, 
subgiant or giant phase. It should be noted that a difference in mass loss rates may exist 
between globular clusters and field stars if only the helium abundance in globular clusters is 
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enhanced, since a higher helium abundance increases the surface temperature of stars and 
therefore shifts the relative position of the instability strip in isochrones. 



VI. SUMMARY 

As we have seen, the possible systematic uncertainties in the helium abundance of glob- 
ular clusters and an instability strip induced mass loss near the turn-off region of globular 
clusters allow an average age of globular clusters as low as ~ 20% below current estimates. 
Therefore the currently quoted age estimate of 14.1±1.5 Gyr [l[ or 15±1.5 @ is poten- 
tially subject to systematic shifts and its uncertainties are model dependent and not fully 
represented by the quoted error. Based on eq. (|5]) and our mass loss calculations, we re- 
calculated the age of five metal-poor clusters using data compiled by Sandage U for four 
cases with and without oxygen-enhancement. Tabulated in table 1, the resultant average 
ages range from 11.3±1 Gyr to 14.6±1.5 Gyr for the four cases, compared with 14.1 Gyr 
from Sandage's estimate assuming a primordial helium abundance, no mass loss, no helium 
diffusion and oxygen-enhancement. In table 1, oxygen-enhancement is approximated by 
models with solar-scaled abundances that satisfy eq. (f|) [p5 ]. In case (iv), the constraint 



from the luminosity function of globular clusters on mass loss holds true with and without 
oxygen-enhancement, and constrains the age of globular clusters to be older than ~ 12 Gyr 
in both cases. 

It is unlikely that all globular clusters will have an enhanced Y of 0.28, given the current 
inferences of the helium abundance in globular clusters. An average Y of 0.26, for example, is 
more compatible with available helium inferences. In all cases, however, after uncertainties in 
the helium abundance, mass loss, helium diffusion and oxygen enhancement are considered, 
and the ~ 10% uncertainty from calibrations of RR Lyrae stars is included, the lowest age 
estimate for globular clusters seems to be about 10 Gyr. 

Clearly, the true accuracy of globular cluster age estimates would benefit from better 
estimates of non-canonical effects in stellar models and better determinations of mass losses 
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in Pop II stars. In inferring the helium abundance, stellar model comparisons with observed 
GC parameters should include non-canonical effects such as helium diffusion, core rotation, 
etc. The existence of an instability strip induced mass loss near the turn-off region can be 
explored by the further observations of the luminosity functions of globular clusters and 
lithium observations in Pop II stars. 

A globular cluster age of 11±1 Gyr is entirely consistent with the age of the beginning 
of nucleosynthesis of heavy elements (i.e., SNe events) measured by radioactive dating that 
gives a lower limit of 10 Gyr [[74,[fq|. It is also consistent with the age of the galactic disk, 
10±2 Gyr, measured from the cooling of white dwarfs ||76|| , provided that the disk collapsed 
from the halo within 1 to 2 Gyr. 

If globular clusters are indeed 11±1 Gyr old, a universe as young as 11 Gyr is then 
possible if globular clusters formed at redshift z > 4. This will significantly relax the 
current constraints on cosmological parameters. For a baryonic universe where 0.01 < Q < 
0.1, the Hubble constant can be as high as 90 km/sec/Mpc; for non-baryonic dark matter 
universes with Q > 0.3, H < 77km/sec/Mpc; if O = 1, H = 60km/sec/Mpc. To turn 
the problem around, if the measurements of the Hubble constant eventually converge, it 
will put constraints on fl in a A = Friedman-Robertson- Walker universe and improve our 
understanding of the age of globular clusters. 
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TABLE I. The age estimates from eq. (5) and our mass loss calculations for five metal-poor 
globular clusters. All ages suffer an additional ~ ±10% uncertainty from the ~ ±0.1 magnitude 
uncertainty in the calibration of RR Lyrae stars. 



Clusters Case (i) a Case (ii) b Case (iii) c Case (iv) d 



[O/Fe]=0 








M68 


12.5 Gyr 


11.3 Gyr 


10.2 Gyr 


M92 


16.2 Gyr 


14.6 Gyr 


13.2 Gyr 


NGC6397 


15.5 Gyr 


14.0 Gyr 


12.7 Gyr 


M15 


14.5 Gyr 


13.1 Gyr 


11.8 Gyr 


M30 


14.3 Gyr 


12.9 Gyr 


11.7 Gyr 


Average 








[O/Fe]=0 


14.6 Gyr 


13.2 Gyr 


11.9 Gyr > 12 Gyr 


[O/Fe]=0.6 


13.9 Gyr 


12.5 Gyr 


11.3 Gyr > 12 Gyr 


a F = 0.24, no mass loss; h Y = 


0.26, no mass loss; C Y 


= 0.28, no mass loss; d an enhanced Y 



and an instability strip induced mass loss according to eq. (15). 
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Figure Captions 

Figure 1. Two isochrones with similar parameters but different mixing lengths. 

Figure 2. Isochrones with and without an enhancement of the helium abundance in GCs 
compared with the fiducial sequence of NGC6397 ([Fe/H]^ —2) (denoted by triangles). The 
Y = 0.24 isochrones are shifted by E(B - V) =0.11 and m — M = 12.35; the Y = 0.28 
isochrone is shifted by E(B — K)=0.14 and m — M = 12.45. 

Figure 3. The propagation of a SN shock front inside a 1O 6 M cloud. 

Figure 4. The calculated turn-off magnitude-age relation for different isochrones. Compar- 
isons can only be made between curves with similar Y, since different Y's yield different 
ZAHB luminosities. 



Figure 5. The calculated luminosity functions for GCs compared with observations ||73|| . iV 
is the number of stars in each 0.2 magnitude bin. The reference point is defined as the point 
on main sequence that is 0.05 mag redder than the turn-off. (a) Y = 0.24; (b) Y = 0.28. 
Mass losses are calculated assuming a mass loss rate of eq. flTSP . The sharpness of the bump 
in subgiant regions is due to the finite grid of stellar models. 

Figure 6. Isochrones with and without a higher helium abundance and mass loss, compared 
with the fiducial sequence of NGC6397 ([Fe/H] ~ —2) (denoted by triangles). The 15 Gyr 
isochrone with Y = 0.24 and no mass loss is shifted by E(B — V) = 0.11 and m — M= 12.35; 
the 12 and 13 Gyr isochrones with Y = 0.28 and a mass loss rate that satisfies eq. (p~5f ) are 
shifted by E(B - V) = 0.14 and m - M = 12.38. 
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